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Abstract 

The whole spectrum of high-energy cosmic ray (HECR) is, very likely, influenced by tachyonic 
neutrinos. Especially, the appearance of two knees can be fitted by the tachyon mass m{v e ) = 
m(vu) — 0.51 eV/c 2 as predicted by a minimal three - flavor model for tachyonic neutrino with one 
parameter 5 = 0.38 eV only. Then the evasion of GZK cutoff could be ascribed to Z°(W ) -burst 
model together with the same mechanism for knees as well as a prediction of left-right polarization 
dependent lifetime asymmetry. A further conclusive experiment might be whether the protons of 
HECR detected on Earth are really right handed polarized? 
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I. INTRODUCTION 



It has been known for years that the observed energy spectrum of primary cosmic rays 
can be well described by an inverse power law in the energy E of charged particles (mainly 
protons) from 10 11 to around 10 20 eV jl, 2] as: 
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where J is the flux in m 2 s x sr 1 . However, the index 7 reveals some abrupt changes as 
follows: 
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Here and E\^ are often called the first and second "knee" while E\Vt}ie "ankle" in the 
cosmic ray spectrum (CRS). In 1966, Greisen, Zatsepin and Kuz'min 3} pointed out that 
nearly no proton with energy higher than Eqzk ~ 5 x 10 19 eV (the GZK cutoff) can reach 
the Earth. The reasons are as follows: A proton will lose its energy via a resonant scattering 
process with a photon in the 2.7 K cosmic microwave background (CMB) radiation: p + 
u 2 7k — ► A* — ► iV + 7r to create a A*, which decays into some baryon N and ir meson. Then 
iV may decay back to a proton again. For every mean free path ~ 6 Mpc of travel, the 
proton loses 20% of its energy on average 4j. A proton produced from a cosmic source (say, 
an active galastic nucleus, AGN) at distance D with an initial E p will on average at Earth 
have only a fraction F ~ (0.8) D ^ 6Mpc of E p . A distance D of 60 Mpc will lead to an order 
of energy degradation (F ~ 10%). So 60 Mpc could be set as a limit within which a proton 
of HECR with energy E p beyond the GZK cutoff can reach the earth. However, no AGN 
sources are known within 60 Mpc of Earth. 

Beginning from 1980's, a number of observational facts and theoretical predictions chal- 
lenged the above claim. Right after the "ankle" energy i% , the cosmic ray spectrum rises 
again and beyond the GZK cutoff [5|]. Although the experimental evidence of these HECRs 
is still a controversial problem today, a possible mechanism called the Z-burst model (ZBM) 



,(2) 



,(3) 
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was proposed first by Weiler {f| and then by Fargion et al. to explain why the GZK 
cutoff could be evaded? In ZBM, it was assumed that an ultra high energy (UHE) neutrino 
v (or antineutrino v) from the remote cosmic space could annihilate with a nonrelativistic 
relic v (or v) in the space within distance D < 60 Mpc to produce a Z° meson: 

u + u-^Z . (3) 

Weiler considered that the neutrino has a nonzero m v so that the resonant process (3) can 
be triggered at a peak energy Er within an interval 5E ~ Tz = 2.495 GeV (Tz, the natural 
width of Z° with mass M z = 91.2 GeV 



3): 



E R = ^ (4) 

If assume m„ ~ leV/c 2 , .Er ~ 4 x 10 21 eV. The resonant cross-section was estimated around 
'. 0~ 23 cm 2 and a hadronic Z-burst would decay to, on average, 30 photons and 2.7 nucleons 

in n 

6j (or 2 protons [7J]) with energy near or above the GZK cutoff. Fargion et al. |7| further 
proposed that the following W + W~ channel reaction chains may be more promising over 
the Z° channel (3): 

v + v^ W + + W~, W ± ^p + X. (5) 

There are two difficulties in the above Z(W)BM. First, the flux of UHE v (or v) may be 
not strong enough (the upper bound on the flux is 5 x 10 _16 m~ 2 s~ 1 sr~ 1 for E u = 10 20 eV). 
Second, the density of relic v (or v) may be too low. The latter difficulty is related to the 
theory of cosmology: Around 1 second after the big bang, when temperature was about 
10 10 K (~ 1 MeV), relativistic neutrinos decoupled and evolved into nonrelativistic relic 
neutrinos today. Now the neutrino temperature is estimated to be T v ~ (^-) 1 / 3 T 7 ~ 1.95 
K (T 7 ~ 2.725 K is the photon temperature of CMB). The mean neutrino (or antineutrino) 
number density for each flavor is (T] Uj ) = (3£(3)/47r 2 )T 3 ~ 54 cm -3 , which seems too low 
for them to be served as target neutrinos in the Z(W)BM. To overcome this difficulty, some 
authors assumed that [9] the relic neutrinos have coalesced into "neutrino cloud", which 
might be a sphere with radius in the range of 10 14 — 10 20 cm (i.e., 7 — 7 x 10 6 au) and has a 



IT] 



believe that 



density in the range of n v ~ 10 10 — 10 16 cm~ 3 . Many other authors, e.g. in 
neutrinos may be a part of so-called hot dark matter (HDM) and may cluster in the halo of 
every galaxy, reaching a density n u = 10 7-9 cm -3 for each flavor. 
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To our understanding, the massive neutrino discussed by previous literatures on Z(W)BM 
is basically a Dirac particle. In this paper we will pose an alternative point of view that 
neutrino is, very likely, a tachyon (whose property will be reviewed in section II). This 
is because the most convincing explanation for two knees in the CRS is focused on the 
tachyonic nature of neutrino (section IV). After a brief discussion on how neutrinos can 
be attached to a neutron star and reach a density about 1.6 x 10 13 cm -3 (for each flavor) 
(section III), we will discuss the kinematics of Z(W)BM in section V rigorously, like that for 
two knees, based on two (momentum and energy) conservation laws. Then in section VI, 
we will claim that the evasion of GZK cutoff might provide the second strong evidence for 
the tachyonic nature of neutrino and could be verified further via possible observations that 
most protons in the HECR should be right-handed polarized. Once this can be tested by 
experiments, the maximum violation of parity conservation law discovered in 1956-1957 will 
be highly emphasized again (see find section VII). Three appendices are added for further 
clarifications. 

II. NEUTRINO AS A TACHYON 

A Dirac particle with rest mass itlq obeys the Dirac equation {h — c — 1): 

i(p = ierVx + m cp, 

(6) 

iX = «crVv? - m x, 

where (p and \ are two-component spinors while cr the Pauli matrices. If one defines 

£ = ^(<P + X), Tj = -y={(f - X) (7) 

Eq. ([6]) can be recast into the Wely representation: 

i£ = 2<xV£ + m r], 

(8) 

it) = — icr'Wri + m £. 

There are two symmetries hidden in the Dirac equation. First, Eq. is invariant under 
the (newly defined) space-time inversion (x — > —x, t — > —t) with: 

<p(-x,-t) -> x(x,t), x(~x,-t) -> <p(x,t). (9) 



Second, Eq. (jSj) is invariant under the pure space inversion (x — > —x, t —*t) 



£(-x, t) -> T](x, t), r](-x, t) £(x, t). 



(10) 



Interesting enough, Eq. ([H]) violates the third symmetry — the (newly defined) pure time 
inversion (x — > x,t — > —t) with 



£{x,-t) ~^T](x,t), T]{X, -t) ^£(x,t), 



(11) 



to maximum. Note that the so-called "time reversal (T)" was introduced to quantum me- 
chanics (QM) since 1932. But actually, the T transformation was a misnomer , it is 
merely a "reversal of motion", including an operation of complex conjugation on the wave 
function (WF) in QM. A "pure time inversion" should be defined clS ( X ^ X j t -t) 



without complex conjugation, it is obviously violated not only in nonrelativistic QM 
but also in the Dirac equation as shown by Eq. (flT|) . 



Ill, 



Now let us compare the Dirac equation with a new equation for neutrino proposed in 



ref. 



12 1 (see also 



13| and Appendix 9B of [ll|): 
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Here the subscripts of fields ^ and r\i [i = e, /i, r) refer to three flavors of neutrinos u e , z/ T , 
they have no mass term originally but are coupled together via one parameter 5. In this 
"minimal three flavors model for tachyonic neutrino" , a neutrino is oscillating among three 
mass-eigenstates because the kinetic relation reads (c = 1): 



E]= V 2 -m% (j = 1,2,3) 



2 a t-2 2 2 r2 

m 1 = 4o , m 2 = m 3 = o . 



(13) 
(14) 



We see that here the neutrino is a tachyon traveling with a group velocity 

(ko dE pc 2 
u = u a = —— = — — = —— > c, (15) 
dk dp E K ' 

whereas its phase velocity 

uj E 

u v = - k =~<c (16) 

Under the basic space-time inversion (Q, Eq. (fl2l) remains invariant after we define 

1 1 

<Pi = ^(& + Vi), x% = -j=(& - m) ( 17 ) 

However, Eq. (12) is no longer invariant under the space inversion (10) (with subscript 
i). This is because the opposite signs in the 5 terms which not only lead to the minus sign in 
Eq. ( Tl3|) . implying the tachyon nature of neutrino, but also endow the neutrino a property 
of maximum violation in parity — there are only left-handed neutrino and right-handed 
antineutrino vf existing in nature (there is no vf and v\ as verified by the experiment [jjj]). 
Interesting enough, these two properties are linked intimately — the permanent polarization 
of a neutrino can be maintained because its velocity u > c whereas an observer's velocity 
v < c. 

More interesting thing is: Eq. (TI21) remains invariant under the pure-time inversion (x — > 
x, t — ► —t) with Eq. (TTTj) (adding subscript i), just relying on the opposite signs in 5 terms. 
Due to the basic space-time invariant Eq. (Q, a pure time inversion will transform a neutrino 
into antineutrino with the same energy but the opposite momenta. We will see its strange 
consequence in the explanation of CRS (section IV and Appendix B). 

Now Eqs. (jl2l) - (116l) are already capable of explaining three experimental facts: 
(a) In the tritium beta decay 3 H — > 3 He + e~ + u e , the neutrino mass square m 2 (defined 
by relation E 2 = p 2 + m 2 ) were reported to be negative in the particle tables of 1996 and 



2000 



li (c= 1): 

m 2 {y e ) = -2.5±3.3eV 2 . (18) 



The pion decay (ti + — > ji + + z/ M ) also shows that [151 ]: 

m 2 (^) = -0.016 ± 0.023 MeV 2 . (19) 

At first sight, one would regard the uncertainty in either (fTBl or ( fl9l) being too big to be 
trusted. However, in the model like Eqs. ( TT2l - (fl6l) . if a neutrino |z/ e ) or |z/ M ) is just created 



6 



from the decay process and staying in three mass-eigenstates, say, 

3 

s > = X)Cik>. (20) 



i=i 



with \vi) having tachyon mass square mf (i = 1, 2, 3), then the neutrino mass square m 2 can 
only take values — 4£ 2 , — 5 2 and — 5 2 with probability ratio |Ci| 2 : |C 2 | 2 : |C 3 | 2 = 2 : 1 : 3 
respectively. If we simply take the average of m 2 , it would be 

(m 2 (z/ e )> = (m 2 (z/ M )> = -2 5 2 ±V2 5 2 , (21) 

where the large "uncertainty" is totally ascribed to the oscillation between — 4S 2 and — 5 2 
even without any errors in measurements. So we suggest that if our experimental physicists 
manage to treat their data in a two-center fitting (with ratio 4 : 1 in (— m 2 ) value and 
statistical weight ratio 1:2), better results could be achieved. 

I — I 

(b) The neutrino oscillation as verified by the Kamiokande group |l7j can be understood 
qualitatively by the solution of Eq. (fT2l : 

\C e (t)\ 2 = l-lsm 2 [±(E 2 -E 1 )t}, 

(22) 

K(t)\ 2 = \CAt)\ 2 = lsm 2 [l(E 2 -E 1 )t}. 
Here |C e (0)| = 1 means that a u e was created inside the sun at an initial time t = 0. 



(c) The accurate data of SNO group 18| showed that the flux of solar neutrinos is likely 
composed of equal amount of three flavors: 

\C e (t)\ 2 -n(t)\ 2 -\Cr(t)\ 2 . (23) 

In our opinion, Eq. (1231) just means that the oscillation among three flavors is not a coherent 
one like that shown by Eq. (I2"2l . but a noncoherent one regardless of the inertial condition. 
This is only possible because the group velocity u g = u > c exceeds the phase velocity 
u g < c{u p Ug = c 2 ), thus destroying the coherence of wave packets and leaving the oscillation 
merely as a mechanism to establish the detailed balancing among three flavors. For further 
discussion, see Appendices B and C. 



III. NEUTRINOS AS A COMPONENT OF THE DARK MATTER 



Before discussing the CRS in following section, let us consider the possibility of neutrinos 
clustering in the cosmic space. According to the theory of cosmology, both neutrinos and 
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antineutrinos are prevailing over the space isotropically. They have an average kinetic energy 
2kT v ~ 3.5 x 1CT 4 eV per particle with a temperature T„ = (^-) 1 / 3 T 7 ~ 2 K. What we try 
to add is two hypotheses: 

(a) They are, most likely, tachyons with a (real and positive) tachyon mass m ~ 0.51 
eV/c 2 (see next section). Under good approximation, the average energy of a tachyon 
neutrino in the cosmos can be evaluated to be 2kT u (rather than ~kT v ). 

(b) The Newton's gravitation law should be complemented according to the invariance 
of mass inversion [m — > —m), which is equivalent to the invariance of space-time inversion 
(at the QM level), showin g th e symmetry between particle and antiparticle — an essence of 



special relativity (SR) [ill. Il9j|. Hence antineutrinos will be repealed by the matter galaxies, 
only neutrinos can be attracted by some celestial body to form a cluster, which might be a 
component of the dark matter. 

For simplicity, consider a celestial body having mass M distributing uniformly in a sphere 
with radius R. Then its gravitational potential reads: 



U{r) 



«M (r>R) 



(24) 



%§(r 2 -3R 2 ). (r<R) 



If a neutrino has kinetic energy 2kT v > \mU(r)\, it will escape from the trap of this celestial 
body. For example, typical star like our sun with M = M & ~ 2 x 10 30 kg, R = R Q ~ 7 x 10 8 
m, f^aiC ^ 1.6 x 1(T 6 eV. A white dwarf has M ~ M , and R ~ iWh ~ 6 x 10 6 m, 



' 2 R 

1.6 x 1CT 4 eV. So either a star or a dwarf cannot capture neutrinos into its trap. 



3GmM ^ i a „ 1 n -4 



2R 



However, a neutron star has M ~ 2M , and R ~ 10 4 m, so ~ 0.23eV > 2kT u . On 

the other hand, though a black hole has R = = fmc 2 ~ |(0.51)eV ~ 0.38 eV, 

no Z° or W created by process (jSJ) or (jHJ) can escape from the black hole. So we will 
not consider it. (If consider a galaxy roughly as a sphere with radius 3 x 10 4 ly and mass 
M ~ 10 n M Q , ~ 2.5 x 10" 7 eV < 2kT u , so we will not consider the neutrino clustering 

in a galaxy either.) 

Let us focus on how many neutrinos N v can be trapped in a neutron star. They can 
spread over the whole volume till a farthest distance r max ~ 2 ^ fc ™ M ~ 6 x 10 6 m > R ~ 10 4 
m. Hence 

Nv = / dN v = / t^4tt / p 2 dp. (25) 



o 



(2vrft) 3 
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Here p max = {[E — mU{r)] 2 + m 2 ^} 1 / 2 ~ {[mf/(r)] 2 + m 2 c 4 } 1//2 since the neutrino energy 
E ~ 0. As ~ 0.212 < 1 for neutron star, we have 

Dividing the integration region into two parts: (0, i?) and (R, r max ), we obtain 

A^, ~ 6 x 10 31 . (26) 

Since neutrinos are trapping basically inside the neutron star, the average number density 
for one flavor is 

n u ~ N„/~irR s ~ 1.6 x 10 13 cm- 3 . (27) 

Many neutron stars are already observed inside our galaxy and they are believed to be one 
kind of the end stage of the star evolution in all galaxies. So the neutrinos trapped inside 
them could be served as targets for UHE antineutrinos incident from the remote space, 
creating Z° and W mesons to escape from the neutron star. However, the total mass of 
these clustering neutrinos seems still quite low. The so-called dark matter in cosmology is 
still a mystery to us. 



IV. EXPLANATION OF TWO KNEES IN THE CRS 



The appearance of two knees in the CRS strongly hint that there is a new reaction channel 



suddenly opened at the threshold energy or E^ respectively as shown in Eq. 



our knowledge, the most attractive mechanism was proposed first by Kostelecky [20] and 



. To 



Ehrlich 



2JJ for the first knee and than refined by Ni for both knees [l3| as follows 



We need only consider a one-dimensional problem, i.e., all processes are constrained 
strictly on a straight line linking the source to Earth (the x axis of S frame) before they 
can be observed by us. consider an UHE proton having velocity v encounters a tachyonic 
neutrino v e with velocity u(> c) in the space (not inside the neutron star because neutrons 
will interact with proton strongly). As u and v are parallel in the S frame, the proton will 
see the neutrino having a velocity u 1 in its rest frame S'\ 

i u — v 
1 — uv/c 2 
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(see Appendix B). Once v exceeds a critical value v cr = c 2 /u(< c), the proton suddenly sees 
the neutrino v e transforming into an antineutrino v e flying toward it with a negative velocity 
u'{< — c). Then a process occurs locally in the S' frame as 

u e + p^n + e + . (29) 

The energy and momentum conservation law for this process are 

m m n m e 
m„ H — = — — H — , (30) 

P V«'7c 2 - 1 v 7 ! " </c 2 Vl-v'l/c 2 

mu' m n v' n | m e v' e ^ 



where m p , m n and m e ( f^, v' n and fg ) are the rest masses (velocities in the S' frame) of 

proton, neutron and positron, respectively. Denoting 

v >. 1 
A = — = tanh C- , 7 4 = , 2 = cosh 0, $7; = sinh C», (» = P, n, e) (32) 
c v 1 - Pj 

with Q being the rapidity of particle i, we have 

m p H -j^=^= = m n cosh + m e cosh (' e , (33) 

m -j== = m„sinhC n + rn, e sinhC e . (34) 

Taking the square of (1331) and (1341) respectively and subtracting them each other, we get: 

cosh(C - O = tt^—H + 2mpm -m 2 n -ml- m 2 ]. (35) 
VSn SeJ 2m n m e [ p ^fff^l 

The constraint cosh£ > 1 leads to 

> (m n + m e ) + m — m . (36) 



Defining 



we find 



2m n m , . 

< 1, (37) 



2 



(m n + m e ) 2 + m 2 — m. 
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> - > 1. (38) 



Rewriting (12B|) as 
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we reveal the condition for the occurrence of an exotic process 



v e + p — > n + e + (40) 

in the S frame as: 

0, (41) 



I > y/l + ^^ + l ^ + l (42) 



or the constraint on the proton energy for the reaction (1401) to occur in the S frame is 



_ m p m p (5 v + 1 _ m p \ p v + E u 

Ep ~ V^W VV " TV Z=k- (43) 

The threshold (minimum) value of _E p occurs at /3„ — > oo or E„ — > 0, it should be identified 
to the energy of the first knee E^, yielding 

rnfi) 1 w \2 2 ii 1.695 x 10 15 Ar , AAS 
E th = 7T~ i( m n + m e 2 + m 2 - m 2 ~ eV, 44 

where m is the tachyon mass of v e in eV/c 2 . Similarly, we consider the threshold energy of 
E p for the occurrence of the process 

+ p -> A + /i + (45) 
in the proton's rest frame (S")as a realization of exotic reaction 

u„ + p -> A + /i + (46) 
in the S 1 frame. Identifying it with the energy at second knee, we find 

,-r(2l 1 r/ n? /9 ?n 3.056 XlO 17 Tr . . 

^ = 2^ [(mA + m ^ + m m ' ] ^ eV ' (47) 

where to/ is the tachyon mass of in eV/c 2 (jtia = 1115.6 MeV/c 2 , m M = 105.7 MeV/c 2 ). 
The observation data of E^ and £^ in Eq. (j5J) show that 

m = m (i/ e ) = 0.54 eV/c 2 to' = to'(i^) = 0.48 eV/c 2 . (48) 

It's a pleasure to see the near equality of m{y e ) ~ m!{y^) as predicted by Eq. ffT2l . which 
gives two eigenvalues of mass for a tachyonic neutrino (either as an outcome of oscillation 
or just created) being 

mi = 25, m 2 (= m 3 ) = 5 (49) 
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with statistical weight ratio 1 : 2, or an average value 

~ ~ Ay/2 

m(v e ) = m(z/ M ) = -5 ± — 5. (50) 

Due to the limited accuracy of present observation, we simply compare the average value of 
flU with dSJJ to find 

5 = 0.38 eV. (51) 

Note that if we consider neutrinos being Dirac particles as in Appendix A, the process 
(I29p or 0451) can occur at much lower energy of E p as long as the velocity u of v is opposite 
to v p and E p becoming higher and higher. So theoretically, in a bath of Dirac antineutrinos, 
there would be no threshold energy for E p to exhibit itself as a knee in the CRS. 

V. Z(W ± )-BURST MODEL BASED ON TACHYONIC NEUTRINOS 

Suppose that an UHE antineutrino v with velocity u annihilates a relic neutrino v with 
velocity u' (in opposite direction with u), creating a Z° meson 

v + v^Z . (52) 



Introducing the notation for tachyons (u > c see 111 ]) 



tanhC = — , sinhC = — =, coshC = — = =, (53) 

u yju 2 /c 2 -l v/u 2 /c 2 -l 

the momentum and energy conservation laws in the S frame read 

mc(cosh C - cosh £') = MVj y/l - V 2 /c 2 , (54) 



mc 2 (sinh ( + sinh (') = Mc 2 /^1 - V 2 /c 2 , (55) 
where M(V) is the rest mass (velocity) of Z°. Combining (|54l) and (|55l) into 

cosh(C + = ^ + 1, (56) 

we expand the left side and denote sinh^ = x, cosh£ = y/l + x 2 , yielding an algebraic 
equation for x: 

x 2 + Bx + C = 0, (57) 
12 



In the solution of (1571) < 1) 



x = smh C ^ -[-(— + 2) smh C ' ± — 2 (l + 2—) cosh C ' 



> 0, 



(59) 



only the plus sign can be used. Then the energy of Z° is 



Ez — E -\- E 



M 2 
2m 2 



2m 2 



(60) 



The lower the energy E 1 of v is, the higher the Ez will be. A highest Ez occurs at E' 
0,p' -> m, yielding (M = 91.2 GeV/c 2 ,m = 0.51 eV/c 2 ) 



M 2 



EZ ' ~ — ~ 8.2 x 10 21 eV 
z 2m 



(61) 



This coincides formally with Eq. (j4j) obtained from subluminal neutrinos (with Dirac mass 
m) as a condition of resonance creation process. 
Let us consider the W /± -burst model 



u + u^ W + + W~. 



(62) 



Based on experience in above discussion, for simplicity, we just consider the situation of 
lowest energy of relic neutrinos E' v — > 0, i.e., u' u — > oo,p' u — > —mc. Using notations for two 
W mesons (1 and 2) like that for n and e + in section IV, we have two conservation laws in 
the 5* frame as 

Mc 2 (cosh Ci + cosh Q 2 ) , (63) 



mc 



^u 2 /c 2 - 1 
mc + 



mil 



Mc (sinh (j + sinh £2) 



(64) 



v/n 2 /c 2 - 1 

where M = 80.4GeV/c 2 is the rest mass of W meson. Combination of fl63l) with fl64l) yields 



m 2 c 4 + 2mc 6 p p = 2M 2 c 4 [l + cosh(0 - C2)]. 



(65) 



When d = (2, i-e., W + and W have the same velocities, the momentum of v reaches its 
minimum, or the threshold energy of incident v for creating W pair is 

2M 2 c 2 



E.-. 



Pu 



2.53 x 10 22 eV. 



rn 



(66) 



We see the minimum energy of W meson being about 1.27 x 10 22 eV, a little higher than 
the maximum energy of Z° meson shown in flBTl) (comparing Fig. 1 of Ref. ?]]). 
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According to jj], being a final product in the chain decay of Z° or W ± mesons, a proton 
p would have energy about (in our notation and estimated value) 



E v v - : 1 x 10 21 eV (for Z* -> 2p + X) 

p 80 80 v r ) 

j-,min 22 

£ P ~ ~ * ~ 3.8 x 10 20 eV (for W ± ^ p + X) (67) 

The estimation of (167D may account for the observed events around 3 x 10 20 eV (above the 
Eqzk ~ 5 x 10 19 eV) and the W + W~ channel into protons may be more promising over the 
Z° channel. 



VI. THREE MECHANISMS RESPONSIBLE FOR THE EVASION OF GZK CUT- 
OFF 

The Z(W)BM provides a possible mechanism for the evasion of GZK cutoff. As explained 
in 7J], those UHE v s may be the pion decay products of initial protons (coming from remote 
AGNs with energy as high as ~ 7.3 x 10 23 eV) via process p + 7 CMB — > p + Nrr, p + 7 CMB — > 
n + Ntc, (N > 2). These UHE protons are impossible to travel from the source to the Earth 
due to interaction with CMB. However, via the complicated decay chain sequence in the 
Z(W)BM, a proton transforms into a u, then Z° or W ± , then a proton again. Thus the 
probability for survival of a proton, though still small (> 1.2 x 10~ 3 ), is at least 120 times 
more favorable than is required for the direct travel from the source to the Earth. 



The second mechanism for the evasion of GZK cutoff was proposed in Refs 2Cj, |21] 



and further stressed in 13] as follows. At the energy above two knees in the CRS, the 
opening of a reaction channel, Eq. (f4"0l) or (146j) . in the S frame provides a new chain of 
decay: p — > n(A) — > p — > n(A) — > p — > ■ ■ ■ . Hence the survival probability of a UHE 
proton is increased considerably. This is because in the decay chain the proton keeps its 



direction within a cone of tiny angle ~ a/1 — (3 2 along its velocity and the n or A nearly 
has no interaction with photons comprising the CMB and the magnetic field in galaxies. 
Interesting enough, a mechanism responsible for deleting proton at energy above two knees 
turns into a mechanism helping to protect UHE protons (above GZK cutoff) from deleting 
by CMB. The previous limit on the distance D < 60 Mpc is now lifted. 

Besides the above two dynamical mechanisms, a third kinematical mechanism was pro- 



posed in Refs 



22 



23 



24 



25j and also stressed in [13( as follows. The lifetimes of n (or 
14 



A) are different for different polarizations: while right-handed n (A) has lifetime r Rh , the 
left-handed one has T Lh 

r r r 

Trh Tlh 1+/?' r jr^w 1 ' 

It is obvious from Fig. [1] that the faster the speed v of n (A) is, the larger the r Rh will 



10 ! r 



3 




0.5 1.0 

Velocity p 



FIG. 1: The lifetime as a function of fermion velocity (3. (a) The lifetime r Rh of right-handed 
polarized fermions. (b) The lifetime r Lh of left-handed polarized fermions. (c) The lifetime r of 
unpolarized fermions. 

be. Hence we expect that the UHE neutron (A) beam in the decay chain will be gradually 
right-handed polarized, so are the UHE protons survived upon arrival on Earth. 

Although further quantitative calculation needs to be done, the combination of above 
three mechanisms, very likely, will considerably enhance the survival probability of UHE 
protons. Hence we are anxious to see the new experimental data further validating the 
evasion of GZK cutoff in the years to come. 
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VII. SUMMARY AND DISCUSSION 



1. It seems very likely that the whole spectrum of cosmic ray with energy ranging from 
10 11 to 10 20 eV, especially the appearance of two knees and the evasion of GZK cutoff, is 
strongly influenced by neutrinos. However, a unified explanation for the whole CRS strongly 
hints that neutrinos are tachyons rather than massive Dirac particles, in conformity with 
other experimental observations in the tritium beta decay, neutrino oscillation and solar 
neutrinos ( see alternative explanation by Wigmans [sol]). 

2. Our simple numerical calculation is based on equal tachyon mass m for neutrinos 
with three different flavors. This equality is verified by Eqs. (I44p . (1471) and (j48j) in fitting 
two knees in the CRS, giving m ~ 0.51 eV/c 2 , which in turn provides strong support to 
the minimal three-flavors model for tachyonic neutrino, Eq. (|12p . containing one parameter 
5 = 0.38 eV only. 

3. Eq. (fT2]) is based on three kinds of symmetries: (a) invariance under the space-time 
inversion (x — > —x, t — > — t); (b) invariance under the pure time inversion (x —>■ x, t —>■ —t); 
(c) maximum violation of pure space inversion (x — > —x,t — > t). The maximum parity 
violation (c) implies that only and i>r exist whereas vr and are strictly forbidden. 
The unique property (b) endows the possibility of a strange transformation of a neutrino in 
the S frame into an antineutrino in another S' frame so that an exotic reaction like (|4T)]) or 
( 146]) can occur for an UHE proton in the cosmic ray. 



4. While the Dirac equation (or reduced Dirac equation, see [26]) shares the same basic 
symmetry (a) with Eq. ffT2"j) . they are just opposite with respect to symmetries (b) and (c). 
Since Dirac discovered his equation in 1928 and wrote it in a four-component covariant form 

as 

(7A + m )^ = 0, (69) 

above three symmetries were hidden and overlooked to some extent for a long time. Even 
we knew the parity is conserved for Dirac equation, it was stated as an invariance under 
the action of an operator P such that ip{x) — > Pip(x) = 74^(— x) ~ where 74 = 

I 

p = I J is a 4 x 4 matrix. In this covariant formalism, it was difficult to understand 

-/ 

the essence of why the parity is conserved in Dirac equation and how it can be violated 
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in some other equation. In our opinion, the two-component noncovariant form for Dirac 
equation, Eq. ([6]), is better in the sense of rendering the above three symmetries accurate 
and explicit as shown by Eq. (jH]) till Eq. (ITT]) . Furthermore, the noncovariant form renders 
the construction of Eq. (FEZ]) possible to be the counterpart of Dirac equation. 

5. The maximum violation of time inversion (t — » —t) for Dirac equation implies that 
we already have a microscopic arrow of time implicitly because we are comprised of Dirac 
particles. So we impose the macroscopic arrow of time explicitly on the description of nature. 
For example, we stay in the S frame and see a proton encountering a left-handed neutrino 
v described by six fields & and rji (i = e, /i, r) 

i 

~ Tji ~ exp [-{p u x - E v t)], (|&| > \r]i\) (70) 

(p u > 0). Once after we jump onto the proton's rest frame S', we will see the v becoming a 
right-handed antineutrino v described by 

^ ~ C c ~ exp [^-(p p x> - E,t% (K c | > |C C |), (71) 

(see Appendix B) with p v < 0. This is because we insist on both t and if flowing forward. 
We never see time flowing backward because the concept of time with arrow is invented by 
human being, not belongs to the particle. In the meantime, we insist on the positive definite 
property for every particle's energy E or its mass m together with the validity of Einstein's 
equation. 

E = mc 2 . (72) 

Therefore, the difference between v and v shown in ( 170]) and ( |7T]) ascribes to the difference 
between i and —i in a strict sense, which is equivalent to the space-time inversion (x — ► 
— x, t — > —t) symmetry, also to the mass inversion (m — > —m) invariance [13|, [19(, showing 
the equal existence of particle versus antiparticle and the essence of special relativity (see 



Appendix C, also [11]). 

6. Physics is essentially an experimental science. Eventually, the validity of the theory for 
tachyonic neutrino can only be tested by experiments. Among others, one further conclusive 
experiment could be the measurement on UHE cosmic ray above the GZK cutoff. If UHE 
protons observed on Earth are really right-handed polarized, they would strongly support 
two things: 

(a) Neutrinos are tachyons with features described by Eq. (TT2l with 5 = 0.34 eV; 

17 



(b) The prediction about left-right polarization dependent lifetime asymmetry, Eq. (1681) . 
is correct. 

These two things in turn imply the maximum parity violation for neutrinos as well as the 
validity of another two symmetries they obey strictly. 
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APPENDIX A: EXPLANATION FOR TWO KNEES AND Z(W)BM BASED ON 
NEUTRINOS BEING DIRAC PARTICLES 

Assume that neutrino is a Dirac particle with rest mass m. Trying to explain the first 
knee in CRS, we can only consider the normal process 

v + p — > n + e + (Al) 

in the S frame. Denoting the velocity of particle % {% = u,p,n,e = e + ) by V{ and 

cosh Q = — 1 sinhQ = — ; V = (A2) 

As in section IV, we consider two conservation laws to get 

2mm p cosh(£ p ± ( p ) = m 2 n + m 2 e — mt — m 2 + 2m n m e cosh(£ n — ^ e ), (A3) 

where "+" or "-" sign refers to v v and v p being "opposite" or "in parallel" . 

(m n + m e ) 2 — ml 1K 

mcosh(Cp ± Q > — = 1-807 x 10 15 eV. (A4) 

2m p 

If consider the v as a relic antineutrino with v D — > 0, i.e., — ^ and m ~ 0.54 eV/c 2 , 

E p = m p c 2 coshC ^ —1.807 x 10 6 eV ~ 3.16 x 10 15 eV. (A5) 

m 

At first sight, the estimated value (A5) looks similar to the energy in Eq. (H4l) for 
tachyonic neutrino. However, they are different due to two distinctions: 
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(a) In (A4) with "+" sign, an increase of (p implies the decrease of ( p . Hence in the 
isotropically distributed relic v bath, (A5) is by no means a threshold value for E p . By 
contrast, Eq. (PH|) is really a minimum energy at E p — > because in Eq. (I4"5I) . E p will 
increase with E v . 

(b) When the reaction (Al) begins to occur for a low energy (i.e., low speed) u, the 
probability can only increase gradually because it is suppressed by the momentum of v, 
Pp — > 0. On the contrary, for a tachyonic v, once the exotic reaction Eq. (j40p is triggered 
at E p = , the reaction probability will suddenly jump out of zero because it is not 
suppressed by the nonzero momentum of u, p v — > m (E v — > 0). 

Next, let us consider the ZBM for neutrino being Dirac particle with rest mass m. Similar 

to calculations from Eq. (|52|) till fl6TT) . but now we denote 

u 1 u I c 

tanhc = — , — = = coshf, — = = sinhc, (A6) 

C y/l - U 2 /C 2 y/l ~U 2 /C 2 



to obtain 



cosh(C + C / ) = ^-l- (A7) 
Ivnr 



where ( and (' refer to the incident v and relic v respectively. Denoting cosh ( = y, we solve 
the equation 

V 2 + B'y + C'^0, 

tf = -(^-2W, C^co^C + (^-^ 
for a root of y = -^K = cosh C and find: 



M 2 

E Z = E, + E U ~— (E' u ±p' v ). ( A8) 



We see theoretically, there is no constraint on Ez- But practically, for the relic neutrino: 
E' v — > m, p'^ — > 0, we have (m ~ 0.51 eV/c 2 ) 

M 2 

^ z ~ = 8.2 x 10 21 eV (A9) 

2m v 7 

coinciding with Eq. (I6T1) for tachyonic neutrino. 

So it seems that for Z(W)BM, the difference between a tachyonic neutrino and a Dirac 

neutrino is not so clear-cut and cannot be discriminated by experimental data available. 

Only by looking at the CRS as a whole, can we claim that the explanation for two knees 

and the evasion of GZK cutoff are stronger evidences for neutrino being a tachyon against 

that for a Dirac neutrino. 
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APPENDIX B: A SOLUTION TO THE TACHYON PARADOX 



The reason why many physicists do not believe in tachyons goes back to a strange puzzle 
involving tachyon motion. See Fig. [2] [ill, Q|- For clarity, we only consider its motion in a 
one dimensional space. 




FIG. 2: A tachyon (P) moving along x axis with velocity u > c. (a) v < t' > 0; (b) v > t' p < 
0. 

A tachyon (P) is moving along the x axis with a velocity u > c in the S frame. We takes 
another S' frame moving relative to S with velocity v. Then if v > c 2 /u, the time coordinate 
of P in the S' frame will become negative: 

t' < (u > c, v > c 2 /u) (Bl) 

which was regarded as the "tachyon traveling backward in time" or "a violation of causality" 



27] . 



In our opinion, the above puzzle can be better displayed in an alternative way. From the 
well known Lorentz transformation (LT), we have the addition law for velocities as: 

= z n (B2) 

1 — uv/cr 

where vl is the velocity of tachyon in the S' frame. As shown in the Fig. [3] [ill . , there is 
a pole at uv/c 2 = 1. For a fixed u, when v increases across the singularity c 2 /u, we will see 
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FIG. 3: Addition of velocities in Lorentz transformation, (a) v! as a function of u for a fixed v; (b) 
u' as a function of v for a fixed u(> c). 

that w' leaps abruptly from oo — > — oo: 

■u' < — c. (m > c 2 /v or v > c 2 /u) (B3) 

However, Eq. flB3j) still remains as a puzzle. According to LT, the momentum p' and energy 
E' of tachyon in the S' frame are related to p and E in the S frame as follows: 

p> = P ~ vE/ ° 2 , E'= E ~ VP , (B4) 

y/l ~V 2 /C 2 ^/1-V 2 /C 2 



with 



ran „ mc 2 



p= , =>0, E = , = > 0. (B5) 
V^/c 2 - 1 v/m 2 /c 2 - 1 



Here m is the tachyon mass of a particle with kinematical relation as 

dE pc 
~dp ~ ~E 



E 2 =p 2 c 2 -m 2 c\ u = d ^ = V 4r>c (B6) 



Combining O with ([B4]) leads to: 



77? 



p = — = — ^^^= (u - v) > mc> 0, (B7) 



i-4J4-i 



2 2 

E' = — , — (c 2 -uv)<0, (u>—orv>—) 

v 2 u 2 1 V U 



Now the puzzle arises: How can a particle have v! < (u > c 2 /v) whereas its p' > ? 



How can it have energy E' < whereas E > ? All of the above puzzles from (1Blj) - (lB7 
comprise the "tachyon paradox". 
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The paradox disappears in a reasonable quantum theory as follows: According to our 
point of view, the tachyon behaves in the S' frame (with v > c 2 /u) just like an antitachyon 
moving at a velocity u' . So its momentum and energy should be measured as: 

p' c = - p > < 0, E' c = —E' > 0. (B8) 

This is because the well known operators in quantum mechanics: 

t=lh i>t (B9) 

are valid only for a particle. For its antiparticle, we should use instead: 

p c = ih—, E c = -ih—, (BIO) 
(where the subscript c refers to an antiparticle) which are just the essence of special relativity 



(SR)[U|. 



Note that, however, the distinction between ( 1B9I) and (1B10I) is merely relative, not abso- 
lute. For example, the energy of positron e + ( which is the antiparticle of electron) in the 
process (1401) is always positive like that of neutron n. But once a neutrino has energy E > 
in the S frame but has E' < in the 5" frame, it behaves just like an antineutrino in the S' 
frame. Then the relations flB8p and flBlOj) must be taken into account. 



APPENDIX C: VIEWING TACHYON AND ANTIPARTICLE VIA ANALYTIC 
CONTINUATION 

We already know that the negative energy state of a particle could be viewed directly as 



its antiparticle's state (see Appendix B, also ll|, |26[]. Moreover, one often views a tachyon 



as the outcome of an analytic continuation from a subluminal particle as follows: 

m c 2 u>c m c 2 m c 2 -im c 2 m s c 2 

E = — > — = = = — = — = E s , (CI J 



where mo — > im s , m§ — > —m 2 ,m s = —irriQ is the tachyon mass being real and positive. 

We will discuss the procedure of analytic continuation more rigorously in mathematics, 
showing that the existence of both tachyon and antiparticle is inferred implicitly but strongly 
by the theory of SR. 
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Consider - = j3 as the real value of a complex variable z and generalize the particle's 
energy into a function of z 

1 



E(z)=m c 2 f(z), f(z) 



(C2) 



a/FWTT)' 

where an extra i has been multiplied. Being a double valued function, f(z) should be single- 
valued on two sheets of a Riemann surface, on which there are two branch points z + = 1 
and Z- = —1, whereas both z = and oo are analytic points, (see, e.g., 
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29|). 




FIG. 4: Two choices of Riemann surface with two sheets (1) and (2) (beneath (1) so invisible) for 
complex variable function f(z). (a) The cut is along the real axis from — oo to —1 and from 1 to 
+oo; (b) The cut is from —I to +1. 

There are two choices of introducing cut to this surface as shown in Fig. H] (a) and (b) . 
Introducing polar coordinates, we write 



z =F 1 = \ z T 1| exp(i#±) = p±exp(i8±) 

i 
2 



/(*) = (p + P-)- 1/2 



exp 



+ 0-) 



(C3) 



The important thing is to fix the rules of how to measure the arguments 9+ and 9- during 
the point z moving along a contour on the Riemann surface as follows: 

(a) When the vector p + or p rotates with z + = +1 or z_ = — 1 as its center, either 
angle 9 + ((f> + ) or #_(0_) is measured from one of the two edges of a cut. If the vector rotates 
clockwise (counter clockwise), the value of angle is negative (positive). 

(b) Once the rotation of a vector is blocked by a cut (on one sheet), we will choose the 
same vector on the opposite edge of this cut and measure its rotation angle beginning from 
zero again. 
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TABLE I: 



'+ 



0_ 



-(<?+ + f-)/2 



/(*) 



implication 



(i) 



P 



Pi 



(2) 



p« = pf 



p(l) _ p(2) 



p(2) _ p (l) 



P 



(1) 



P 



(2) 



7T 



— 7T 



-7T 



7T 



7T 



7T 



— 7T 

— 7T 

— 7T 
7T 




— 7T 
7T 



7T 



-7T 



-l 



l 



Pi 



Pr 



Pr 
Pl 



(c) The analytic continuation means that we can follow the motion of variable z (or its 
vector), through the gap and enter the second sheet (2) from the first one (1) (or vice versa). 
During this process, the value of argument change should be accumulated in counting. 

We start from the point P[ , located at the real axis {(3 > 0) of sheet (1) inside the 
unit circle on Fig. 4 (a) and bring the z point along a contour until it reaches P 2 located 
at the upper edge of cut. Because two sheets of Riemann surface are cross-connected via 
this cut, the function value of P 2 i s equal to that of P 3 [on the sheet (2)] which is just 
beneath the Pg at the opposite (lower) edge of cut on sheet (1). So the motion of z point 
can be continued along the doted line (clockwise) until it reaches the point P{ which is 
just beneath the P x . 

By the rules mentioned above, we obtain the Table HI and HT1 The last column shows what 
kind of particle (with spin |) may be inferred as a consequence of analytic continuation. 
While p and p imply a subluminal particle (including the proton, of course) and its antipar- 
ticle, v and v imply a superluminal particle and its antiparticle (because the neutrino may 
be the only tachyon we know today). The subscript L or R refers to the particle's helicity. 
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TABLE II: 



-(<£+ + 0_)/2 



/(*) 



implication 



(i) 



Of } = Q? ] 



— IT 

-2tt 
-3vr 

-47T 




-2vr 









TT 






Pl 



TT 



3tt 
2 



-1 



-2vr 



Pl 



TT 



i 

-i 



To study the tachyon more clearly, we perform a transformation z — > w = - with 



/(*) - /(-) 



9{w) 



-i)w 



W y/(w-l)(w + l) 

Now the particle's energy becomes a function on the Riemann surface of g(w): 
E(w 



m c 2 g(w) = m c 2 r (r + r_) x ' 2 x exp< i 



a ( 

2 



! 2 



(C4) 



(C5) 



where w = roe ta ,w =F 1 = r±exp(i0±). There are also two choices to introduce the cut as 
shown in Fig. 5 (a) and (b). The outcome of analytic continuations are listed on Table HTT1 
andlVl 

In Fig. |4] (a), all points describing a subluminal particle p are located within the interval 
(— 1, 1) on the first sheet (1). Its corresponding state of antiparticle p is just located at the 
same point on the second sheet (2). In Fig. [5] (b), p and p are located at opposite edges 
of the cut. This may be of no surprise. The strange thing is: In the outside region with 
\z\ > 1 occupied by tachyons v and u, they are not separated into different parts but living 
in mixing. For example, even P 2 an d -P5 belong to the same upper half sheet (1) so can 
be reached each other by an analytic continuation, they describe v and v respectively. 
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TABLE III: 



w 


0+ 




a 




- + 0_)-§ g(w) 


implication 


Q?> = Qt ] 





— IT 










1 


Pr 




7T 


— IT 







7T 

2 








TT 


— IT 


TT 




7T 

2 


i 


v R 


Q ( i ] = 





TT 







— 7T 


-1 


Pl 


Q? = Q? 


TT 





TT 







1 


Pl 


Q? = Q? } 


— TT 





— TT 




— TT 


-1 
Vl"/3 2 


Pr 


Qf 


— TT 


— TT 







TT 

2 




"r 


Q? ] 


— IT 


— TT 


— TT 




7T 

2 


—i 
V/3 2 -l 




TABLE IV: 


w 


0+ 


4>- 


a 


a — t;(4>. 


f +</>-)- 




implication 



1 

i 
i 

—i 

1 

\/l Z /3 2 ~ 
-1 

-1 



5(1) 



-7T 



Sd) = B(2) 



7T 



5(1) _ 5(2) 



-2tt 



5(1) _ 5(2) 

r Z ~~ ^2 



-2tt 



— TT 



5(1) 



7T 



7T 



5(2) 



-3tt 



5(2) 



-7T 



-7T 



26 




FIG. 5: Two choices of Riemann surface for complex variable function g(w). (a) The cut is made 
in favor of tachyons inside the unit circle: \w\ < 1; (b) The cut is made in favor of subluminal 
particles outside the unit circle : \w\ > 1. Points pj~ (Q^) are counterparts of pP\Q^ ) in Fig. [31 



To understand it, we turn to Fig. [5] (a) where tachyons occupy the whole real axis inside 
the analytic region \ w\ < 1. To our surprise, while Qjp and Q% describe u and u respectively, 
they are living so close and even can transform into each other via the analytic point w = 



Some of our readers may wonder how far can we go by pure mathematics. Let us look 
back at Fig. 4 (a), what is the difference between two states of a particle p (with spin 1/2) 
with (3 > and (3 < 0? It's nothing but its helicity. This is because during the analytic 
continuation only the particle's velocity is changing smoothly without reversal of its spin 
orientation in space. Once when a particle is at rest {(3 = 0), there is no distinction between 
its helicity being left-handed or right-handed at all! 

The unique feature of u L and v R is: They have opposite helicities and are particle and 
antiparticle as well. How can they transform into each other when (3 — > oo? Let us look at 
their wave function as follows: 



i\) v ~ exp 



h 



(p u x - E v t) 



—i 



(p p x - Ept) 



(C6) 



They are really approaching each other because E v ~ E v ~ e — >• and p u ~ — p p , meaning 
that a v L is capable of transforming abruptly into a v R moving at the opposite direction 
when their energies approach to zero. 

We even cannot prevent ourselves from guessing that such kind of process might occur 
somewhere in the universe. 

Actually, similar possibility is already explored to explain two knees in the CRS as dis- 
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cussed in this paper. Unlike the UHE proton, we can easily chase a v L with velocity u — > oo 
and see it transforming into a v R flying toward us once our velocity v exceeds a small critical 
value v cr = c 2 /u [see Fig. [3] (b) and Table IV] . 
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TABLE V: Comparison between two kinds of fermions — Dirac particle and neutrino 

Dirac particle neutrino 
Symmetry under the space- invariant invariant 

time inversion, Eq. 



Symmetry under the pure 
space inversion, Eq. (flQl) 



Symmetry under the pure 
time inversion, Eq. (j lip 
particle velocity 
(group velocity) 



Oscillation among various 
flavor eigenstates 



dynamical equation 
of particle 



invariant 
(parity conservation) 

noninvariant (maximum 
violation — a time arrow) 
subluminal 

(u < c) 



impossible 
(strictly forbidden) 

particle with each flavor 
(like e, r) has its own Dirac 
equation with different rest mass 



noninvariant 
(maximum parity violation) 
[see Eq. (JT2D] 
invariant 



superluminal 
(u > c) 
subluminal (u p < c) 
possible (incoherent oscillation 
leading to distribution with equal 
probabilities among three flavors) 
neutrinos with three flavors 
share a same Eq. (]12p with 
only one coupling constant 5 



phase velocity superluminal (u„ > c) 
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